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We obtain the standard quadrature-phase positive operator-valued measure (POVM) for homo- 
dyne detection directly and rigorously from the POVM for photon counting without directly employ- 
ing the mean field approximation for the local oscillator. In addition we obtain correction terms for 
the quadrature-phase POVM that are applicable for relatively weak local oscillator field strengths 
and typical signal states. 
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I. INTRODUCTION 



With the advent of squeezed states of light Q , a full quantum description of optical homodyne detection Q, |^ H, Q 
assumed importance as homodyne detection (HD) yields phase-dependent measurements of the light field. Whereas 
photodetectors acquire phase-insensitive information about photon statistics 0, 0| , homodyne detection mixes the 
signal field with a coherent local oscillator (LO) to yield photon statistics on the output fields that depend on the 
phase Lp of the LO. By varying this phase Lp, phase-dependent properties of the signal state p can be inferred. The 
phase-sensitive measurement with respect to the in~phase quadrature x or its canonically conjugate out-of-phase 
quadrature p, or some in-between quadrature x^p = Pi^_^/2 = xcostp + psmip, is necessary to observe the nonclassical 
properties of squeezed light. Phase-sensitive measurement has developed beyond measuring specific quadrature-phase 
statistics to acquiring information for many values of and reconstructing the density matrix p for the signal field. 
This technique, known as optical homodyne tomography 0, illustrates another important application of homodyne 
detection. Homodyne detection has evolved into a key tool of quantum optics with applications including squeezed 
light detection, optical homodyne tomography and continuous variable quantum teleportation |^ |^ [lol | . 

The homodyne detection scheme discussed above involves mixing the signal field with a LO field at a beamsplitter 
(BS), and the two output fields are subjected to photodetection, as shown in Fig. The measured photodetection 
statistics are analyzed to infer the quadrature-phase statistics. Only in the limit of infinite LO field strengths can 
the measurement be said to correspond to quadrature-phase measurements, and, of course, this limit is in principle 
unattainable. However, a good approximation to quadrature-phase measurements is attained. In the most useful 
variant, a 50/50 BS is used, and the difference between the photocounts at the two output ports is used to infer 
the quadrature-phase statistics. This is known as balanced homodyne detection (BHD) and has the advantage of 
automatically canceling the photon number sum at the two input ports from the detected output fields. 

The description of homodyne detection begins with photodetection of the output fields and then, to validate the 
approximations normally applied in homodyne detection of quadrature-phase POVM, must show that the resul- 
tant two-mode photon statistic reduces in some way to the quadrature-phase distribution, for the signal field p. 
This connection between photon statistics to quadrature-phase, or joint quadrature-phase, measurements has been 
established via calculations involving quasi-probability distributions or characteristic functions (moment-generating 
functions) for the electromagnetic field and allowing the local oscillator strength to become infinitely large. Yuen and 
Shapiro introduced the characteristic function approach in their seminal quantum theory of HD j3] , and Walker em- 
ploys Wigner functions in his analysis of HD ■ Braunstein uses the positive P-representation in the description 
of the photon counting statistics, and he emphasizes the quantum nature of the LO as he investigates "the effects of 
a finite-amplitude fully-quantum-mechanical local oscillator" . Banaszek and Wodkiewicz [T^ calculate moments of 
operationally defined quadrature operators, with an emphasis on finite photodetection efficiency, but in contrast to 
our approach, employ the mean field approximation to the LO from the outset. 

These studies undoubtedly establish the connection between the exact photodetection statistics and the approx- 
imate quadrature-phase HD. However, modern applications of homodyne detection, for example to quantum infor- 
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FIG. 1: Balanced homodyne detection scheme: (a) the input state I?/;) is mixed with a LO in coherent state |q), and photon 
counting occurs at the two output ports; (b) photon numbers j ± m' are shown entering the two input ports and j ± m are 
counted at the output. 



mation applications such as continuous-variable quantum teleportation, requires an operational quantum theoretic 
approach [iJI; Banaszek and Wodkiewicz advocate the operational approach, but here we avoid the mean field ap- 
proximation and thereby include correction terms for the POVM corresponding to HD. The operational approach is 
important in the context that a measurement may be applied for some purpose other than characterizing the state /5; 
paradoxically, in continuous-variable quantum teleportation |^, ^3 j the sender mixes the field described by density 
operator p with one component of a two-mode squeezed vacuum state |l5j in such a way that the sender cannot 
know, even in principle, what the density operator p was that is being subjected to this measurement. For such 
applications, a rigorous approach to homodyne measurement, which demonstrates that the POVM for photodctection 
reduces to the POVM for quadrature-phase or joint quadrature-phase measurements is necessary. Here we establish 
this connection between actual and convenient POVMs by directly calculating the photon counting probabilities using 
two different approaches: (i) working in the Fock basis for the Hilbert space of the signal and LO modes, we employ 
asymptotic expressions for SU(2) Wigner functions that are the BS matrix elements in the Fock basis; (ii) working 
in the over-complete basis of coherent states and taking advantage of the simple transformation of coherent states at 
the BS, we employ the Glauber-Sudarshan P-function. 



II. BALANCED HOMODYNE DETECTION SCHEME 



A balanced homodyne detection scheme is depicted in Fig.^a). The (generally mixed) signal state p to be measured 
is coherently mixed at the BS with a LO assumed to be in a coherent state (in the optical domain a coherent state 
with an absolute adjustable phase has not been achieved, but the coherent state approach leads to correct measured 
results provided that the signal field and LO field are derived from the same source 0, The photon number 

difference from the two BS output ports is measured. The photon number sum can also be measured but usually is 
not. However, in our analysis we include the treatment of both the difference and the sum as this is a more complete 
description than considering the difference alone. We will denote the photon number difference by 2m S Z and the 
sum by 2j € Z. 

The Hilbert space of two modes of electromagnetic field has the basis {|ni) (g) \n2)} of joint eigenstates of the 
photon number operators hi = a\ai and fi2 = 0,20,2- Denoting j — (rii + ri2)/2, m ~ {rii — n2)/2, we will use the 
notation \ jm) = \ni) <Si \n2)- Thus, the state \jm) is the number states with photon numbers j ± m at modes 1 and 2, 
respectively. The value of j can be any non- negative half- integer and m can get values — j, ~j + 1, . . . , j for a given j. 

A. Beam splitter transformation 

The beam splitter action on a two-mode state of electromagnetic field is given by the SU(2) transformation [Hill 

m 

5(^1,(^2,(^3) = e-*'^i-^^e-'^^-^«e-'^^^% (1) 
where the SU(2) generators Jy, Jz are expressed in the Schwinger boson representation as 

Jy ^ -]^{a\a,2~ a\a,i), = ^(ajai - 4^2) (2) 
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An input state \ jm) is transformed under the BS action as 

m' 

where m' in the sum runs from — j to j with unit steps and dl^,^{ip2) — {jm'\e~'^'^'^'^^\jm) are the SU(2) Wigner 
functions [20| . 

On the other hand, coherent states are transformed in a very simple way on BS. If the initial two-mode coherent 
state is lai) Cg) |a2), where 



a,) =e ^1")' (4) 



?i=0 

then the BS output state is again a two-mode coherent state with amplitudes a'l, o^^'- 

B{ipi,^2, ® 1^2) = \a[) (g) 14) (5) 

with a[ = e"''^i/2(Q;i cos<p2e"''^='/^ - 012 sin^sac''^^^^) and a'2 = e''^i/2(Q;i sin(y92C"''^-^/^ + 02 cos (^se'"^-^^^). This sim- 
ple transformation is a key reason for the usefulness of the Glauber-Sudarshan P function in describing homodyne 
detection. 

For the rest of the paper, we will consider BHD with no phase factors, so we set (pi = Q,(p2 = 7i'/2, 1^33 = 0. Ideally, 
the LO is prepared in the pure coherent state with amplitude a, and is directed into port 1 of BS. The unknown 
signal field described by the density operator p enters the second input port. The total state of the two modes before 
entering BS is then 

Pin = \a){a\ (g) p (6) 

The beam splitter transforms the input state into 

Pout = BpinB^ (7) 

The probability of detecting j + m and j — m photons at the two BS outputs is then 

PI = Tr(pout \jm) Om|) = {jm\BpinB^ \jm) (8) 
The probability can be expressed as 

PL = Tr(/5outEi), (9) 

where the POVM = \jm){jm\ satisfies the completeness condition J^mj^m = ^ and positivity condition 

Tr(pi,iE:^J > 0. The importance of the photon number difference measurement and its relation to phase measurements 
has been emphasized for many years, including in early work on phase operators in two mode systems pl| . If the 
total photon sum 2j is not measured in BHD and only the difference 2m is observed, the appropriate POVM is 

00 

Em = \jm){jm\ 

j = \m\,\m\ + l,... 
00 

= ^ \j + m)i{j + m\(g)\j ~7Ji)2{j -ml, (10) 

J — I m I , I m I + 1 , . . . 

where the subscripts refer to the output ports of BS. However, we consider here the more valuable case when both m 
and j are measured. 



III. ASYMPTOTIC SU(2) WIGNER FUNCTION APPROACH 

In this section we derive the photon counting probability in the strong LO limit using the asymptotic formula 
for SU(2) Wigner functions. Let p be the density operator describing the signal state and let the coherent amplitude 
of the LO be a = —A with A real and positive. We now prove the following Theorem: 
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Theorem 1 For A very large (in the limit A—fOo), which means a very strong LO, the photon counting probability 
P-!^ is given by 

^-(2j-A^f/2A^ 

PL - — — (^l/5|2;>, (11) 

where \x) is the eigenstate of the quadrature operator x ^ (a + a))/\/2 with the eigenvalue x = mj \f]. 

Proof: We assume a pure signal state first, p = with j-i/;) = X^^o''^" !"■)■ 1"^^ photon counting probability 

is given by the square of the magnitude of the probability amplitude for j ± m photons emerging from the 
first/second interferometer output port, 

Ml = Om|S(0, V2,0)| - = e-^'/^ ^ ^„ ^77=^ '^^„.,-„('^/2). (12) 

V(2J-")! 



ri=0 



For A — > oo, the probability distribution of the total photon number 2j is dominated by the Poissonian distribution 
of the photon number in the LO, so 2j is sharply peaked at . Further, the photon number difference 2m at the 
BS output is much less than 2j and also n <JC 2j holds for any photon number n for which -!/;„ is non-negligible. This 
enables us to use several approximations. First, the fraction in Eq. H12|) can be approximated via the Stirling formula 
and the Taylor expansion and by neglecting terms of order n? / j, nj j and higher. These approximations yield 



V(2j-«)! W 

The condition n <C 2j justifies the following asymptotic expression for rf^j-ra(^/2) that holds for n <^ j and is 
central to the calculation: 



<,-„(V2) = (-l)"r'/'^n(v<?arcsin"' 



J 

« i-irr'/'^u^im/^). (14) 

Here Un{x) = {x\n) denotes the n*^ Hermite Gaussian, that is, the x-representation of the number state The 
approximation -y/Jarcsin(TO/j) sa rn/y/j is valid for \m\ <^ j. Substituting Eqs. and (|14|1 into Eq. H12|l and 
approximating j by /2 in the denominator, we obtain 

Ml^ = ^- ^Vnw„(m/yj). (15) 

^''^ n=0 

The following identity for the inner product of the state and the eigenstate holds due to completeness of the 
Fock basis: 

oo oo 

{x\il]) = ^ {x\n){n\ip) = ^ u„(a;) V'n- (16) 

n=0 n=0 

Eq. Hlt)|) also holds if the summation over n goes only to 2j instead of infinity because 2j ^ n for all n for which ijjn 
differs from zero significantly. Then Eq. (|15|l becomes 

ML = ^ (xH) (17) 

with the eigenvalue x = m/y/j. Eq. Ullfl is now obtained directly by squaring the magnitude of ML for the pure 
signal state. The extension to mixed states is straightforward and follows from linearity of quantum mechanics. □ 
In the case of a general phase of LO when the amplitude is a = — Ae"^, Eq. Hll() turns into 

PL = ^/^A vi^\p\^)v^ (18) 

where \x)^ is the eigenstate of the rotated quadrature — xcosip +psmip with the eigenvalue x — m/^/j. 
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From Theorem n] we can now get the POVM defined by Tr(/3ii:JjJ = P-!^ and corresponding to BHD in the strong 
LO hmit: 

n^- ^- \X){X\. (19) 

Eqs. and (|19() show that in the limit of strong LO, homodyne detection performs the POVM given by the 
projection |a:;)(a;| to the x-eigenstate. This fact has been known; however, here it has been shown for the first time 
by a direct calculation. However, our result does not provide any correction terms. We will obtain these in the next 
section by employing the Glauber-Sudarshan P function. Before doing so, let us discuss a few aspects of the result 

m- 

First, the Gaussian factor in Eq. reflects the fact that the Poissonian distribution of the photon number for 
the LO (whence the majority of the total 2j photons come) converges asymptotically to the Gaussian distribution 
P(2j) = (27r)-i/2^-i exp[-(2j - ^2^/2^2]^ 

Second, one may wonder if the probability distribution is properly normalized. Indeed, it is easy to check that 

00 

E E = ^ (20) 

2j=0 m.= -j,-j + l,...j 

by changing the double sum into an integral and using the normalization of the state p, 

/oo 
{x\p\x)dx^l, (21) 
-00 

and replacing x — m/y/j by a; = \/2m/A, which can be done for a strong LO. 

Third, if the total photon sum 2j is not measured in the homodyne detection scheme, then the probability distri- 
bution for the photon number difference 2m is 

00 

i^™ = E ^™ ^ (^1^1^) (22) 

j = \m\ 

(in the sum j runs from \m\ to infinity via unit steps and the eigenvalue x is again V2m/A). The factor 1/V2A 
in Eq. H22|) is connected with the Jacobian V2/A of the map m ^ x = \^m/A and the fact that m changes in 
half-integer steps. 

IV. GLAUBER-SUDARSHAN P-FUNCTION APPROACH 

The method using the asymptotic formula; for SU(2) Wigner functions from the previous section gave us the 
asymptotic expression for the photon counting probability P,^. However, it is difficult to obtain the correction terms 
due to the amplitude of the LO being finite because of absence of correction terms in Eq. H14|) . This problem can be 
overcome by using the Glauber-Sudarshan coherent-state representation, which we do in the following. 

We represent the signal state p by the Glauber-Sudarshan P-function |^.|2^|2^ 



^ P(/3)|/3)(/3|d2/3. (23) 

The BS input state is then 

Pin= |a>(a|® / P(/3)|/3)(/3|d^/3, (24) 



and the BS output state is 

a — f3\ / a — (3 



Pout = / P(/3) 



V2 /i\ ^/2 



a + P \ / a + P 
V2 /aX \/2 



d^p. (25) 
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Using Eq. the probability P^^ is evaluated as 

{j + to) ! ( j — to) ! 



P(/3) e-l'^l' \a - |q, + /JpO-™) d^/?. (26) 



We again assume that the LO amplitude is a = —A. Generalization to arbitrary a is straightforward and discussed 
later. 

To evaluate the integral in Eq. H26|l , we use the following identity, definitions and lemma that is proved in Appendix 

A: 

Identity 1 For |a;| < 1 



(1 + x)" = exp[nln(l + x)] = exp 



. k=l 



Definition 1 A pure z-regular state is a state that can be expressed in the Fock basis as 

oo „ 

|^)=AA5]^|n) (27) 



n—0 ^ 



with the complex coefficients Cn satisfying |c„| < 1, z €! R"'" and Af a constant. In other words, it is a state whose 
Fock basis coefficients fall off at least as fast as those of a coherent state \z). 

Definition 2 A mixed z-regular state is a finite mixture of pure z-regular states, that is, a state corresponding to 
density operator 



Y.Pr\A){H (28) 



with n finite, Pi > and all \tpi) being z-regular. 

Examples of z-regular states include (i) a coherent state I7) with I7I < z, (ii) superposition or mixture of several 
such coherent states, (iii) superposition of such a coherent state with a number state, (iv) superpositions or mixtures 
of several number states. However, they do not include squeezed or thermal states. 

Lemma 1 The Glauber- Sudarshan P function P{P) of a z-regular state is identically equal to zero for \(3\ > z. 

(See Appendix A for the proof.) 

We assume that the signal state is z-regular for some z < A. Then P(/3) = for \f3\/A > 1, and we can employ 
Identitynin evaluating the powers \a + and \a — /3|2(j-™) in Eq. as follows: 

\ k=l fc=l / 

Using Eq. (|29() . the integral in Eq. (|26|l can be expressed as 



where the exponent 

X = -^^{/3 +(/?)} + -E2fc^ ^1 ^Efc I2I ■ (31) 

k=2 k=2 

To evaluate the integral H3U|) . we will use the following lemma. 
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Lemma 2 Let p = J ^'(7)|7) (7I d^7 be the Glauber-Sudarshan representation of the density operator p. Then for 



^(7)7™(7*)"e 
(For the proof see Appendix 0). 



'(7+7*)-^]= = V^Tr {p ia^r\x){x\ a™) 



(32) 



Theorem 2 For a z-regular state p and the LO coherent amplitude —A with A > z, 




^/2? 



\/2? 



A 



A 



1 



+ («t) 



t\2fe-l 



2fc- 1 



^2fc-l 



^ 1 a^fc + (at)2'.^ 

fe=2 



^2/c 



(33) 



The ordering symbol : : that involves the projection operator |a;)(a:;| should be understood as 

: \x){x\a'-ia^Y ■.= {a^Y\x)(x\a\ (34) 

that is, all creation operators go to the left of the projector |x)(a;| and all annihilation operators go to the right of it. 

Proof: The theorem is proved by a straightforward calculation applying Lemma El to Eq. (|30|l and substituting the 
result into Eg. □ 

The form of in Eq. (|33() produces the POVM for homodyne detection of a z-regular state: 



W = 





^/27 



A 



2.7 - A"^ 
2^2 



{a2 + (at)2} + 2m^ 



1 



~2fc-l 



t^2fc-l 



k=2 



2fc- 1 



^ 1 a2fc + (at)2fe 

k=2 



A2k 



(35) 



such that Tr{pn^„} = holds. 

Eq. (|35|) is the key result of our calculation. It shows that the POVM for homodyne detection of a z-regular state 
(with z < ^) is given by the normally ordered product of the projector multiplied by an exponential of powers 

of creation and annihilation operators. We will discuss this result in the following sections. 

We still need to mention the case of a general phase of the LO when a = —Ae^'^. The operators a and in Eqs. H33|l 
and P5|l then have to be replaced by = e"^a, a}^ = e~"^a^, respectively, and \x) has to be replaced by \x)^, the 
eigenstate of the rotated quadrature x^p ~ x cos ip + psin f = {a^ + aJ,)/\/2 with the eigenvalue x = \/2m/ A. 

A. Limit ^ — > 00 

We begin discussing the result H33|) by considering the limit of strong LO, that is, the limit A ^ 00 for a given 
signal state p. This will give us the asymptotic expression for the photon counting probability P/„ corresponding to 
an ideal homodyne detection. 

For large A, the total photon number distribution is dominated by the Poissonian LO distribution, so 2j is peaked 
at A^ and has the variance of A^ . Hence the expression (2j — A^)/2A^ in the exponent of Eq. (|33|) is negligible. At the 
same time, in the sums in the exponent the factors A~^ go to zero for A 00. Thus the trace in Eq. H33() becomes 
simply Tr(p|x)(a;|) — {x\fi\x). The factor in front of the trace can be approximated using the Stirling formula for the 
factorials and neglecting terms of order {2j — A^)/2A^ and mjA?!"^. Then the probability P^^ becomes 



PL 



RA^ 



{x\p\x), 



(36) 
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which rephcates the result Hll|l from the Sec. lIIII The only difference is that in Eq. Hll() the eigenvalue was x = m/\/j 
while here we have x = ^f^mjA. However, this difference is not important as j is sharply peaked about for a 

strong LO as has been mentioned. 



B. Infinite series and its convergence 

For a finite amplitude of LO, one can expand the exponential function in Eq. H33|) using the usual Taylor series. 
This gives an expansion of the the photon counting probability into the following scries: 

V^2-2^ e-^'A4^ e2™VA^ J 2.7 - 

PL - — + — y^ApV) - [^""l" ^'^^ + '\ApW) k>] 

+ |^[(x-|a3p|x) + (x|p(at)3|:,)] + ...|. (37) 

The terms in the series are arranged such as to contain increasing powers of creation and annihilation operators. To 
determine for which states this series converges is a task that we have not been able to solve in general. We believe, 
though, that the following conjecture is valid: 

Conjecture 1 The series in Eq. 15*71 ) converges for all z-regular states with z < A. 

Surprisingly enough, however, it turns out that the question of convergence does not really matter for practical 
purposes as we will see in the following section. 

In addition, also the factor in front of the parentheses in Eqs. or (|37|) can be expanded into a series using 
the Stirling formula for the factorials and Taylor expansion around the point m/j = and (2j — A'^)/2A^ = 0. The 
leading term of the series for this factor is the same fraction as in Eq. H36|l and reflects the Gaussian limit of the 
Poissonian distribution for the LO photon number. We do not write the other terms explicitly. 



C. Truncation in Fock basis 



We explore the properties of the series (|37|l for density operators truncated in the Fock basis. Such density operators 
can be expressed as 



N 

/5= ^ Prnn\m){n\ (38) 
m.n—0 

for some finite N. 

Theorem 3 For truncated signal states p, the series \'Jl\j is finite (i.e., it contains only a finite number of non-zero 
terms). Therefore it converges and expresses the exact photon counting probability P^^. 

Proof: Consider a term in the series in Eq. I|37|) that contains more than 2N field operators (i.e., annihilation and 
creation operators). Then it contains more than N creation and/or more than iV annihilation operators. As all the 
annihilation operators are to the left from the density operator p and all creation operators are to the right of it, every 
such term turns into zero because of the truncation (|38|l of p. Further, it follows from the expansion of an exponential 
in Eq. H33|) that in the series in Eq. H37|l the number of terms with less than k field operators is finite for every k. 
Hence, the number of nonzero terms in the series in Eq. H37|) is finite, which wc wanted to prove. □ 

The fact that the series converges for truncated states is very useful as it can be employed for states for which the 
series does not converge. The reason is the following. Consider a general state p = n=o'°™" ^"^^ fo'^ ^ 

given cutoff € N define the corresponding truncated state p' with matrix elements pj„„ satisfying 

Pmn tor m < N,n < N ^^^^ 
otherwise 

This definition ensures the proper normalization of p'. Now, the cutoff number N can be chosen arbitrarily large, 
so that the truncated state p' mimics the state p arbitrarily close. Then also the photon counting probabilities P^' 
corresponding to the state p' can be brought arbitrarily close to the probabilities P^^ for all pairs of j, to, for which 
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is non-negligible. This enables us to employ Eq. H37() for calculating P/^^ with an arbitrary precision also for states, 
for which the series (|37|1 does not even converge. 

Another question concerns the practical usefulness of this truncation procedure. To see an example when it is not 
useful, consider the signal state as a coherent state with an amplitude /3, |/?| 3> A, and its truncation for a very large 
N (say iV ^ |/3p). In this situation, the series H37() diverges while after the truncation it becomes finite and so it 
converges. A closer inspection of Eq. H37|l also shows that the initial subsequent terms grow very quickly for both the 
original and truncated states. Therefore wc would need very many of them to calculate the probabilities P,jj using 
the truncation procedure and Eq. H37|) . which would not be very practical and it would be much simpler to calculate 
directly. This can be expected as the signal field is not weaker than the LO field. 

On the other hand, in many situations our result is very useful. Our calculations were motivated by trying to 
show that homodyne detection measures the field quadrature, and to find correction terms. This happens for large 
amplitudes of LO when the term (a;|p|x) in the photon counting probability P^ is the largest and dominant one. In 
such situations the truncation works very well and the series (|37|l gives good correction terms for balanced homodyne 
detection as can be seen in the following section. 

We should also note that the convergence of the series (|37(l is not directly related to the behavior of the initial terms. 
It can happen (e.g. for a weak thermal state or a weakly squeezed vacuum state) that the initial subsequent terms 
decrease quickly but after some time, they start to grow and the scries diverges. At the same time, for weak signal 
states (compared to the LO) these first terms provide an increasingly good approximation to the photon counting 
probability P^ as can be seen in the next section with numerical simulations. The situation is thus similar to the 
one in perturbation theory: even though a perturbation series diverges, its several (or many) initial terms may give 
a good approximation. 



D. What is a strong local oscillator? 

We would like to address the question now of when the LO is strong enough so that BHD really performs the 
projective measurement of the quadrature phase of the signal field. It can be roughly said that it is in situations for 
which the first term {x\p\x) in the brackets in Eq. (|37|) dominates over the remaining ones. Let us focus at the second 
and third terms, 



[{x\ci'p\x) + {x\p{a^f\x)], (40) 



2A2 



and try to estimate their magnitude compared to (x|p|x) . First, the distribution of the LO photon number is Poissonian 
with both mean and variance equal to A^. Therefore, if we assume that the LO contains many more photons than 
the signal state, the quantity {2j — A'^)/2A'^ is of order of 1/A. Of course, 2j can be an arbitrary integer, but if it is 
not close enough to A^, the probability P^ becomes negligible. In this sense we mean that {2j — A'^)/2A'^ is of order 
of 1/A 

To estimate {x\a'^p\x) + c.c, we will consider two different types of signal states - a coherent state and a number 
state. The discussion for a general state would be very difficult, and we think that coherent and number states are 
good representatives that can help us understand the general behavior of the series in Eq. H37() . 

For a coherent state |/3) for which p = |/3)(/3|, 



-2 



(a; a 



p\x) + {x\p{a^f\x) = 2Re{/3^} {x\p\x). (41) 



This means that the term (|40|) in the series (|37ll is of order of Ke{P^}/A compared to the first term {x\p\x). We see 
that if the mean photon number in the signal state is much less than the magnitude of the LO amplitude, the leading 
term is dominant. 

For the signal field in a number state \n) we have 



(xla^plx) + {x\p{a^f\x) = 2y/n{n- l)Re{{x\n){x\n - 2)*}. (42) 

The magnitude of the inner product {x\n — 2)* can be considered roughly the same as that of {x\n)* for our purpose. 
As n{n ~ 1) is close to n for n > 1, we arrive at a similar result as for the coherent state: the second and third 
terms become unimportant if A is much larger than the photon number in the signal state. 

The analysis of the magnitude of other terms in Eq. 137|l would be similar. The result is that if ^ ^ n, where n 
means the average photon number in the signal state, the subsequent terms decrease quickly and homodyne detection 
indeed measures the field quadrature phase. It should be noted that it is not enough if the mean number of photons 
n in the signal state is much less than the number of photons in the LO; in fact, the correct condition is that the 
square of n must be much smaller than the number of photons in the LO. 
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This condition has a clear physical interpretation. As the photon number in the coherent state | — ^) has a Poissonian 
distribution of width A, the condition of a strong LO can be formulated such that the mean photon number in the 
signal must be much less than the fluctuation of the photon number in LO. Now suppose for a moment that the 
opposite would hold. Then from the knowledge of the total photon number 2j we could access some information 
about the photon number in the signal state. However, the photon number operator does not commute with the 
quadrature x, so this would necessarily disturb the measurement of x. On the other hand, if the strong LO condition 
is satisfied, then we do not know how many of the 2j photons come from the signal and how many come from the 
LO; thus, the different possibilities can interfere and the distribution of x is not affected. 



E. Numerical simulations 



In this section we show some numerical simulations of our results. For a given pure signal state {tp) and a given 
photon number sum 2j we compare the exact photon counting probability calculated with the help of Eq. H12(l 
with the series (|37|) truncated at different points. The purpose of such a simulation is to show that taking increasing 
number of terms in the series H37|) gives an increasingly better approximation to the exact probability . 

The LO amplitude was chosen to be a = —A = —20 which means that the mean photon number of the LO field 
is 400. The value of j in the individual plots was chosen randomly from the Poissonian distribution of LO photon 
number. It has turned out during the simulations that changing j inside the interval for which the probability P^ is 
non-negligible does not affect the behavior of the series significantly. As the signal states we have chosen a coherent 
state with amplitude 2, a squeezed vacuum state exp[r(a^ — a^^)/2] |0) with r = 1.5 and a number state |6). The 
results of the simulations are shown in Fig. El The exact probabilities P^ are shown in black, and the results of 
truncation of the series (|37|l keeping terms with (i) zero number of field operators, 

PI''' = ^ „ (^IpN), (43) 

(j + my. (j - my. 

are shown in green color, (ii) maximum of two field operators 

= ^ [m^) - -^Ma'p\^) + (-lp(at)^lx)]} (44) 

are shown in blue color, and (iii) maximum of four field operators are shown in red (we do not write Pm^' explicitly). 

The simulations show that with increasing number of terms in the series 1)3 7f) . a better approximation to the exact 
photon counting probability is achieved. 



V. CONCLUSION 



We have analyzed balanced homodyne detection in terms of the POVM for photon counting by directly calculating 
the photon counting probability. We employed two different approaches. First, using asymptotic expressions for SU(2) 
Wigner functions allowed us to establish the non-trivial connection between the discrete variables j, m corresponding 
to photon numbers being detected and the continuous quadrature phase variable x^. In the strong LO limit, we 
have shown that homodyne detection indeed performs the projective measurements corresponding to POVM |a;)(a:|, 
where \x) is the eigenstate of quadrature phase operator. Second, employing the Glauber-Sudarshan P-function, we 
extended the result obtained by the first approach. For a very large amplitude of the LO, the result was the same, 
and for finite amplitudes we obtained additional correction terms. Even though the series we got does not converge 
in general, it can be used for determining the correction terms via truncation of the signal state in the Fock basis. 
We have determined the strong LO condition for coherent and number signal states - the square of the mean photon 
number in the signal state must be much smaller than the mean photon number in the LO. We have also performed 
numerical simulations that confirm the validity of the quadrature-phase POVM and the correction terms for a LO 
that is not strong for typical signal states. Therefore, in addition to obtaining the quadrature-phase POVM rigorously 
from the photon counting POVM, we have an expansion that yields correction terms for the POVM that works well 
for typical signal states in quantum optics. 

In this paper we have considered a perfect HD scheme with ideal detectors, LO and BS and 100% mode-matching. 
In practice, all these elements are subject to imperfections, which disturbs the measurement. For example, the LO 
from a realistic laser has an amplitude distribution Plo(q^) broader than the delta-function. This would convolute 
the probability P^(q;), where we now write the dependence on a explicitly, with Plo(Q!)- If this distribution is 
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Gaussian, then the strong LO measurement would correspond to a Gaussian spread of the quadrature measurement 
with the imprecision corresponding to the degree of LO amphtude fluctuation. Lossy beamsphtters and inefficient 
photodetectors would add vacuum noise that would result in Gaussian spread of quadrature measurement, similar to 
the effect discussed above for the LO amplitude spread. Finally, for a multi-mode field with the LO mode-matching 
condition satisfied, the detection efficiency can incorporate the mismatch between the beam and detector modes. If, 
on the other hand, the signal and LO modes are mismatched, HD efficiency declines, and beats between different 
frequency modes arise. 

Our operational approach to HD ignores the realistic effects described above, but the theory is readily generalized 
to accommodate these effects by including inefficiencies and multimode description. Moreover, if multimode fields 
and beats are desirable, heterodyne detection replaces homodyne detection (for which signal and LO are frequency 
matched) ; an operational formulation of heterodyne detection without mean field approximation is a topic of further 
research. 
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APPENDIX A: PROPERTIES OF THE P REPRESENTATION FOR z-REGULAR STATES 

We first prove Lemma ^ for pure z-regular states and then generalize to mixed states. The density operator 
associated with a normalized pure state is 

p=\^l^){^\, (Al) 

and is represented by the Glauber-Sudarshan P representation according to Eq. (|23|l . The trace of p is unity due to 
the normalization of the state j^/') and hence the integral of the P-function over the complex plane is equal to unity: 



P(/3) d'p = J P(/3)Tr 1/3) (/3| d^/3 = Tr p = 1. (A2) 
Now, for a positive number r we define a non-unitary operator 

S{r) = e-^ exp (A3) 
and, for a normahzed z-regular state lip) [see Eq. (|77|) ]. we consider the state 

°° / n 

|V'')=^WK'>=AA^^|n), (A4) 
where the coefficients are related to the coefficients c„ by 



m=0 



Clearly |c^| < 1, so the state {ip') is also z regular but generally not normalized. To normalize it, we introduce the 
inverse norm A/"' = {tp'l'ip')-^^^ so that the state IV'') jv = .U'lip') = ('0'l'0')~^^^l'0') is normalized. The density operator 
p' of \4'')n can be expressed via the density operator p as 

p' = cxp p exp (—\ (A6) 
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and the P function corresponding to \tp') n is hence 
p'(/3) =AA''e-''- cxp 



exp 



r/3* 
z 



P{(3) = Af' exp 



2r 



Rc/3 



P(/3). 



The integral of P' over the complex plane is unity as the state \iIj')n is normalized: 

P'(/3)d2/3= 1 



Decomposing /3 to real and imaginary parts /3 = /3i + i/32, and using Eq. IjATp . we can write Eq. 
integral 



(A7) 

(A8) 
as a double 



— oo 
oo 



where we have denoted 



d/3i exp 
exp 2r 

G(/3i) 



2.(^-1 



d/32P(/3i,/32) 



G(/3i)d/3i 



P(/3i,/?2)d/32 



2n 



(A9) 



(AlO) 



(All) 



(A12) 



We see that the integral (jA12p is bound by a fixed number AA^e^^ , no matter how large r we choose. The only way 
to satisfy this is if G{(3i) = V/3i > z. Specifically if f{x) : R — > R is a function and wc know that 

(A13) 



(A14) 



The inner product {ilj'\ip') can be bound as follows [see Eq. (|A4|l ]: 

n=0 ' n=0 

From Eqs. (|X9)l and I|A11|I . 

G(/3i)d/3i < AA^e"'. 



poo 

/ exp 




' —OO 





/(.t) e'^^ dx < c 

where c > is fixed and a is arbitrary positive, then necessarily f{x) = for all x > 0. Thus, we obtain 



/oo 
P(/3i,/?2)d/32 = for /3i>z, 
-oo 



which means that the integral of P over any vertical line in the complex plane that is farther than z from the origin 
is zero. 

Now the whole construction can be repeated with another state 



Wu>) = e ''exp 



whose P function is 



p;(/3) = (^;i^;)-iexp 



2r 



(A15) 



(A16) 



Using the same argument, we arrive at the fact that the integral of P{P) over any line whose normal has the angle ip 
with the real axis and whose distance from the origin is larger than z (eg the line / in Fig. [2Jl is zero. Now, as cp can 
be arbitrary, this means that the integral over all lines not intersecting the circle with radius z is zero. Then it follows 
by the tomographic argument that the P-function must be zero outside the circle, which is what we wanted to prove. 

The generalization of the claim to mixed z regular states is straightforward as the P-function of a mixed state is 
the weighed sum of the P-functions of the pure states in the mixture. 
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Im p 



Rel3 



FIG. 2: The integral of the P function of the state |?/)) over the line I shown in the picture vanishes, and so it does for all other 
lines that do not intersect the circle |/?| = z. Therefore the P function itself vanishes outside the circle, i.e., for |/3| > z. 



Here the fact that |(a;|7)p ~ tt ^/^exp[— — \/2Rc{7})^] was used. 
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APPENDIX B: PROOF OF LEMMA |21 



A direct calculation yields 
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FIG. 3: The simulation for (a) coherent state I7) with 7 = 2 for j = 190, (b) squeezed state with squeezing parameter r — 1.5 
for j = 219.5, and (c) number state |6) for j = 183.5. The exact probabilities are shown in black, and the truncated ones are 
shown in green, blue and red, respectively, according to the increasing number of terms in Eq. I|37|l taken into account. The 
red curves are so close to the black ones in (b) and (c) that they almost cover them in the plots. 



